/*=========================================================================
The Impact of Denying a Wanted Abortion on Women and Children
Authors: Juliana Londoño-Vélez and Estefanía Saravia

Creation date: January 9th, 2024
--------------------------------------------------------------------------
Table 2 - Baseline Characteristics of Women Filing Abortion Rights Claims
=========================================================================*/

use "$data/Final_Medellin_SisbenIII.dta", replace

	
global outcomes Edad_Encuesta Edad_Tutela Embarazo_adolescente  Nivel_0 Primaria Secundaria Media Superior puntaje_float Estrato_0_1 thogar Any_Children Numero_Hijos  Soltera Casada Separada Viuda Cabecera
	
mat a=J(19,9,.)

matrix rownames a =  $outcomes "N"


local i=1
foreach var in  Edad_Encuesta{
	* (1)
	sum `var' 
	mat a[`i',1]=`r(mean)'
	mat a[19,1]=r(N)
		
	* (2)
	sum `var' if tutela==1
	mat a[`i',2]=`r(mean)'
	mat a[19,2]=r(N)
	
	* (3)
	sum `var' if Caso_Aborto==1
	mat a[`i',3]=`r(mean)'
	mat a[19,3]=r(N)
		
	* (4)
	sum `var' if Caso_Aborto==1 & Niega==1
	mat a[`i',4]=`r(mean)'	
	mat a[19,4]=r(N)
	
	* (5)
	sum `var' if Caso_Aborto==1 & Niega==0
	mat a[`i',5]=`r(mean)'	
	mat a[19,5]=r(N)

	* (6)
	gen comp_`var'=`var'*Niega if Caso_Aborto==1
	ivreghdfe comp_`var' (Niega=Mujer_Ponente) if Caso_Aborto==1, by(oficina#time) cluster(id_juez)
	mat a[`i',6]=_b[Niega]	
	
	* (7)
	ttest `var', by(tutela)
	mat a[`i',7]=r(p)
	
	* (8)
	ttest `var', by(Caso_Aborto)
	mat a[`i',8]=r(p)
	
	* (9)
	ttest `var' if Caso_Aborto==1, by(Niega)
	mat a[`i',9]=r(p)

	local i=`i'+1
	}

	
		local i=2
foreach var in  Edad_Tutela Embarazo_adolescente  {

	* (3)
	sum `var' if Caso_Aborto==1
	mat a[`i',3]=`r(mean)'

	* (4)
	sum `var' if Caso_Aborto==1 & Niega==1
	mat a[`i',4]=`r(mean)'	

	* (5)
	sum `var' if Caso_Aborto==1 & Niega==0
	mat a[`i',5]=`r(mean)'	

	* (6)
	gen comp_`var'=`var'*Niega if Caso_Aborto==1
	ivreghdfe comp_`var' (Niega=Mujer_Ponente) if Caso_Aborto==1, by(oficina#time) cluster(id_juez)
	mat a[`i',6]=_b[Niega]	

	* (9)
	ttest `var' if Caso_Aborto==1, by(Niega)
	mat a[`i',9]=r(p)
	
	local i=`i'+1	
	}

	
local i=4
foreach var in  Nivel_0 Primaria Secundaria Media Superior puntaje_float Estrato_0_1 thogar  Any_Children Numero_Hijos Soltera Casada Separada Viuda Cabecera {
	
	* (1)
	sum `var' 
	mat a[`i',1]=`r(mean)'
	mat a[19,1]=r(N)

	* (2)
	sum `var' if tutela==1
	mat a[`i',2]=`r(mean)'
	mat a[19,2]=r(N)

	* (3)
	sum `var' if Caso_Aborto==1
	mat a[`i',3]=`r(mean)'
	mat a[19,3]=r(N)

	* (4)
	sum `var' if Caso_Aborto==1 & Niega==1
	mat a[`i',4]=`r(mean)'	
	mat a[19,4]=r(N)
	
	* (5)	
	sum `var' if Caso_Aborto==1 & Niega==0
	mat a[`i',5]=`r(mean)'	
	mat a[19,5]=r(N)
	
	* (6)	
	gen comp_`var'=`var'*Niega if Caso_Aborto==1
	ivreghdfe comp_`var' (Niega=Mujer_Ponente) if Caso_Aborto==1, by(oficina#time) cluster(id_juez)
	mat a[`i',6]=_b[Niega]	
	
	* (7)
	ttest `var', by(tutela)
	mat a[`i',7]=r(p)

	* (8)
	ttest `var', by(Caso_Aborto)
	mat a[`i',8]=r(p)
	
	* (9)
	ttest `var' if Caso_Aborto==1, by(Niega)
	mat a[`i',9]=r(p)

	local i=`i'+1
	}


	svmat2 a, rnames(outcome)
	
	drop if outcome==""
	
	rename a1 mean_sisben
	rename a2 mean_anytutela
	rename a3 mean_abortion
	rename a4 mean_abortion_denied
	rename a5 mean_abortion_nondenied
	rename a6 mean_compliers
	rename a7 p_value_1_2
	rename a8 p_value_1_3
	rename a9 p_value_1_5
	
	
	keep outcome mean_* p_value*
	
	order outcome mean_* p_value*
	
	export excel using "$output/Tables paper matrices replication.xlsx", sheet("Table 2", modify) 
	